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Abstract 

A broad extension of periodic system into the sector of antimatter could be possible sometimes in 
a remote future. We expect that anti-alpha spontaneous emission from an antimatter nucleus will 
have the same Q-value and half-life as alpha emission from the corresponding mirror nucleus. This 
is the consequence of the invariance of binding energy as well as of the surface and Coulomb energy 
when passing from matter to antimatter nuclei with the the same mass number and the same 
atomic number. The Q-values and half-lives of all measured up to now 27 cluster radioactivities 
are given together with Q-values and half-lives of the most important competitor — a decay. The 
lightest anti-alpha emitter, 8 Be, will have a very short half-life of about 81.9 • 10 -18 s. 
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I. INTRODUCTION 


In 1928 Dirac predicted the existence of negative energy states of electrons [lj when he 
developed his famous relativistic wave equation for massive fermions. The antimatter char¬ 
acter of these states became clear in 1933 after discovery of the positron (the anticlectron) 
in cosmic radiation by Anderson 


Individual anti-particles are produced by particle accelerators and in some types of ra¬ 
dioactive decay. Antiprotons ( p ) 3] were observed in 1955 by Segre and Chamberlain. The 
antineutron was discovered in proton-proton collisions at the Bevatron (Lawrence Berkeley 
National Laboratory) by Cork et al. in 1956 ( 4 ] . Antiprotons are produced at Fermilab for 
collider physics operations in the Tevatron. Other accelerators with complex projects for 
antimatter physics are the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory, LHC at CERN, and in the future Facility for Antiproton and Ion Research 
(FAIR)’s high-energy storage ring in Darmstadt f|. 


Until now it was established in all experiments that every antiparticle has the same mass 
with its particle counterpart; they differ essentially by the sign of electric charge. Also every 
antinucleus has the same mass or binding energy as its mirror nucleus j(3]. 


Anti-atoms are difficult to produce; the simplest one — the antihydrogen ( H) was pro¬ 
duced, cooled and confined [7] for about 1000 s js-ll]. At the beginning the Low Energy 


Antiproton Ring (LEAR) at CERN was used. This device decelerated the antiprotons and 
stored them in a ring. The antimatter helium-4 nucleus, 4 He, or anti-cq consists of two 
antiprotons and two antineutrons (baryon number B = —4) [l2j. This is the heaviest ob¬ 
served antinucleus to date. It seems that the next one, antilithium, has an extremely low 
production rate. 


It will be a long way to produce a rich diversity of more complex antinuclei justifying 
a broad extension of periodic system into the sector of antimatter and strangeness 13]. 
Nevertheless in this work we try to understand whether their decay modes by anti-a and 
anti-cluster spontaneous emission would differ from a decay and cluster radioactivity 
of corresponding mirror nuclei. 
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II. POTENTIAL BARRIERS 


Let us assume that a binary decay mode (e.g. anti-alpha decay, anti-cluster decay or 
spontaneous fission) of a parent anti-nucleus, A Z , leads to an emitted anti-cluster, Ae Z e) 
and a daughter anti-nucleus, Ad Z d - 

A Z -A Ad Z d + A *Z e (1) 

with conservation of baryon numbers. Alternatively the subscript d may be denoted with 
1 and e with 2. By definition, the number of antiprotons of an antinucleus is equal with 
the number of protons of the corresponding nucleus. The same is true for the number of 
antineutrons and of neutrons. Consequently, there is a good reason to assume that every 
anti-cluster, Ae Z e , will have the same binding energy as the cluster Ae Z e , and similarly the 
binding energy of the daughter anti-nucleus, Ad Z d , will be identical with that of the daughter, 
Ad Z d , and the binding energy of the parent anti-nucleus, A Z, is identical with that of the 
parent A Z. The released energy 


Q = [M-(M d + M e )]c 2 


( 2 ) 


can be calculated using the last evaluation of experimental data for atomic masses 17]. In 
this eq. c is the light velocity, ill, M d , M e are the masses of parent, daughter and emitted 
nucleus. 

In general the ratio of Z/A ^ Z d /A d ^ Z e /A e meaning that the three partners have 
different charge densities. One can take into consideration the difference in charge densities 
El by assuming uniformity in each of the two fragments. In this way the nuclear volume 
V = V\ + V 2 is divided in two parts, each of them being homogeneously charged with a 
density 

' r e 1/1 (3) 

p 2 e, r G V 2 

During the decay process from one parent to two fragments there is a potential bar¬ 
rier which determines the metastability of any anti-nucleus. It is penetrated by quantum 
mechanical tunnelling as was shown by Gamow in 1928 for alpha decay of nuclides 19]. 

For cylindrical symmetry the simplest parametrization of the shape during this process, 
with only one deformation parameter (the volume and the radius of the emitted fragment 


Pe{ r) = 


3 






are conserved), is that of two intersected spheres assumed in the two-center shell model 


20 ]. The radius of the initial spherical anti-nucleus is Ro = r^A 1 ^ and the radii of the two 
fragments are R e = r 0 A^ 3 and Rd = r^A 1 ^. Within Myers-Swiatecki’s liquid drop model 


(LDM) 
(Y+EM) 


the radius constant r 0 = 1.2249 fm and in the Yukawa-plus-Exponential model 


22] r 0 = 1.16 fm. During the overlapping stage the separation distance of the 


two fragments increases from an initial value R t = Ro — R e to the touching point value 
R t = R e + Rd- II is convenient to use the deformation parameter £ = (R — R i )/(R t — Ri ) 
equal to unity at the touching point R = R t . 


We apply the macroscopic-microscopic method [23] to calculate the deformation energy, 
Ed e f, according to which a small shell and pairing correction SE is added to the macroscopic 
phenomenological model deformation energy obtained by summing the surface and Coulomb 
energy due to the strong and electrostatic forces: 


J def 


= (E s - E ° s ) + (E c - E° c ) 


( 4 ) 


where E° s = a 2 o7l 2 / 3 = a s (l — k s I 2 )A 2 ^ 3 and Eq = 3e 2 Z 2 /(5r 0 7l 1//3 correspond to the spherical 
parent with a s = 17.9439 MeV, I = (N — Z)/A and k s = 1.7826 within LDM. 

The proton levels and neutron levels of a single particle shell model, e.g. two center 
shell model 24], allowing to calculate { 23 ] the shell and pairing correction, SE, are different 
because protons are electrically charged. In the same way for antinuclei the antiproton levels 
should be different from antineutron levels but the antiproton levels would be identical with 
proton levels and antineutron levels identical with neutron levels. 


A. Strong interaction 


During the deformation from R = Ri to R = R t the strong interaction is responsible for 
the surface energy. The strong force acts between antinucleons in the same manner it acts 
between nucleons; the electric charge doesn’t play any role. For a number of antinucleons 
equal to that of nucleons it will have the same effect. The deformation dependent term is 
obtained by division with E° s \ 


_ _ a 2 i a 22 

d s — — — - n s 1 H- n s2 

Eg a 2 0 a 2 o 


( 5 ) 
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with (Z 21 7 ^ 022 7 ^ O 20 taking into account the difference in charge densities, ihn and B s2 are 
proportional with surface areas of the fragments: 


B s i — 


d 2 r x 
2 " i-i 


1 

4 


V' 

dx 


1/2 


dx 




' 2,1 w' 2 ’ 
y hi 


11/2 


dx 


( 6 ) 


(7) 


where d = (z" — z')/2Rq is the length of the deformed antinucleus divided by the diameter 
of the spherical shape and x c is the position of separation plane between fragments with - 1 , 
+1 intercepts on the symmetry axis (surface equation y = y(x) or yi = y(x')). 


B. Coulomb interaction 


We can see that not only the surface energy but also the Coulomb energy is invariau 


when passing from matter to antimatter because in the following general relationship 
the charge density appears as a product of p e (r)p e (ri): 

E = 1 f f p e (r)p e (r 1 )d 3 rd 3 r 1 

c 2 J Vn 7 | r — iq | 

See also the expression of Eq above. 

For fragments with different charge densities by dividing with E^ we obtain 


25 1 


( 8 ) 


B c = H = (—\b c1 + + CW B a 

E c \P0eJ Pte \POeJ 


(9) 


explicitly showing the electrostatic self-energies and the interaction of two fragments. For 
binary systems with different charge densities and axially-symmetric shapes, we got 


B c i = b c 7 dx J dx'F(x, x') 

/ Xc rl 

dx / dx 1 F(x, x') 

~1 1 J Xn 


( 10 ) 

(11) 

( 12 ) 


B c2 = b c dx dx'F(x , x') 

J X c * 7 X c 

where b c = 5 d 5 / 87 r and d, x c were defined in the previous subsection. The integrand is given 


by 


F(x,x') = {yyi- 
+K 


K -2D 

2 (y 2 + yl) - (* - 

- x 1 ) 2 + ^(x — x') 

( 

3 

[ 

If* 

0 x — x' dy 2 

V 2 dx _ 

2 | x-x'dyl 
Vl 2 dx' 

1 


dy 1 _ dy 2 ' 
dx' dx 


>} a P 


-1 


(13) 
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where D = (K — K')/k 2 ] K and K' are the complete elliptic integrals of the first and second 
kind, respectively: 

K(k) = I (1 — k 2 sm 2 t)- l l 2 dt ; K'(k) = j (1 — k 2 su\ 2 t) 1 ^ 2 dt (14) 

Jo Jo 

and a 2 = (y + yi) 2 + (x — x') 2 , k 2 = Ayyi/a 2 . The elliptic integrals are calculated by using the 
Chebyshev polynomial approximation. For x = x' the function F(x,x') is not determined. 
In this case, after removing the indetermination, we get F(x,x’) = 4y 3 /3. 


C. Examples 


In figures [0 and [2] (top panel) we present two examples of potential barriers calculated 
within LDM and Y+EM for spontaneous emission of 34 Si from 242 Cm and 14 C radioactivity 
of 250 C/, respectively. In both cases it is clear that within Y+EM the strong interaction 
continues to act, as a proximity force even for separated fragments, R > R t , as long as the 
the tip separation distance remains small enough; the interaction energy is maximum at 
certain distance R m > Rt- For spherical fragments there is an analytical relationships of 
interaction term: 


Ev 12 — —4 


Ri 


r o 


\J + 21+22 


Ri, 


9i92 4 


R 


— 92 fl — 9lf2 


exp(— R/a) 


R/a 


g k = — cosh — - sinh — ; f k = — 


Ri 


Rl 


sinh ( — 


(15) 


(16) 


where a = 0.68 fin is the diffusivity parameter and a 2 = a s (l — ft/ 2 ), a s = 21.18466 MeV, 
ft = 2.345. 

The contribution of surface, E s , and Coulomb energy, Ec, to the LDM potential barrier 
is plotted at the bottom of figures |T] and [21 The potential barrier height is the result of 
adding an increasing with separation distance surface energy up to the touching point with 
a decreasing electrostatic energy up to infinity. 


III. HALF-LIVES 


The experimental data on halflives against cluster radioactivity (261, [27|, T c , and a decay, 
T a , are given in Table [0 together with Q-values, updated using the mass tables published 


in 2012 


17|. Up to now there was not observed any odd-odd cluster emitter. 
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TABLE I. (Rvalues in MeV and decimal logarithm of the half-lives in seconds for the most probable 
CR and aD of cluster emitters experimentally observed. 


Parent 

Emitted 

Qc 

log l0 T™r(s) 

log w T c ASAF (s) 

Qa 

log w T™r(s) 

22i Fr 

14 C 

31.291 

14.52 

14.27 

6.458 

2.55 

221 Ra 

14 C 

32.395 

13.39 

13.74 

6.881 

1.90 

222 Ra 

14 C 

33.049 

11.01 

11.15 

6.679 

1.58 

223 Ra 

14 C 

31.828 

15.19 

14.72 

5.979 

5.99 

224 Ra 

14 C 

30.534 

15.86 

15.93 

5.789 

5.50 

226 Ra 

14 C 

28.196 

21.19 

20.98 

4.870 

10.70 

223 Ac 

14 C 

33.064 

12.96 

12.68 

6.783 

2.48 

225 Ac 

14 C 

30.476 

17.28 

17.69 

5.935 

6.23 

228 Th 

20 0 

44.724 

20.72 

21.72 

5.520 

7.78 

231 Pa 

23 F 

51.860 

26.02 

25.52 

5.149 

11.47 

230-g 

22 Ne 

61.386 

19.57 

20.12 

5.992 

6.26 

230 Th 

24 Ne 

57.760 

24.61 

24.86 

4.770 

12.38 

231 Pa 

24 Ne 

60.409 

23.23 

23.01 

5.149 

11.47 

232-JJ 

24 Ne 

62.309 

20.42 

20.37 

5.413 

9.34 

233 F 

24 Ne 

60.485 

24.84 

24.97 

4.909 

12.78 

234 U 

24 Ne 

58.824 

25.92 

25.72 

4.857 

13.04 

235-jj 

24 Ne 

57.362 

27.42 

29.97 

4.678 

16.57 

233-g 

25 Ne 

60.727 

24.84 

25.48 

4.909 

12.78 

235pj 

25 Ne 

57.706 

27.42 

30.38 

4.678 

16.57 

234 U 

26 Ne 

59.415 

25.92 

26.59 

4.857 

13.04 

234 U 

28 Mg 

74.109 

25.14 

25.34 

4.857 

13.04 

236 Pu 

28 Mg 

79.668 

21.52 

20.55 

5.867 

7.95 

238 Pu 

28 Mg 

75.909 

25.70 

25.60 

5.593 

9.44 

236-g 

30 Mg 

72.274 

27.58 

29.54 

4.573 

14.99 

238 Pu 

30 Mg 

76.795 

25.70 

25.86 

5.593 

9.44 

238 Pu 

32 Si 

91.186 

25.27 

25.33 

5.593 

9.44 

242 Cm 

34 Si 

96.509 

23.15 

22.77 

6.216 

7.15 
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It is clear that T c » T a , hence cluster radioactivity of nuclei with atomic numbers 
Z = 87 — 96 is a rare phenomenon in a huge background of a particles. The measurements 


are m goo c 


model 


28 


agreement with predictions within analytical superasymmetric fission (ASAF) 
Surprisingly, for some super heavy nuclei we found jit], 30] comparable half- 


lives or even shorter T c < T a . 

We expect that the same Q-values and half-lives will be observed in the future for anti- 
cluster decay and anti-alpha decay of antimatter nuclei. Perhaps the easiest way to observe 
the decay modes of antimatter nuclei would be to produce the lightest a emiter, s Be, which 
will be split in two 4 He or two a nuclei with a half-life of about 81.9 • 1CU 18 s= 81.9 as - 
the same with that of 8 Be —> a + a 31]. 
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FIG. 1. (Color online) TOP: potential barrier for emission of 34 5i from 2i2 Cm calculated within 
LDM (red) Y+EM (black). BOTTOM: two main terms of the LDM barrier: surface energy (dashed 
line cyan) and Coulomb energy (dotted line blue). 
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FIG. 2. (Color online) TOP: potential barrier for emission of 14 C from 250 Cf calculated within 
LDM (red) Y+EM (black). BOTTOM: two main terms of the LDM barrier: surface energy (dashed 
line cyan) and Coulomb energy (dotted line blue). 
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